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Abstract 

The Relativistic Mean Field (RMF) theory is applied to the analysis of ground- 
state properties of Ne, Na, CI and Ar isotopes. In particular, we study the recently 
established proton skin in Ar isotopes and neutron skin in Na isotopes as a function 
of the difference between the proton and the neutron separation energy. We use the 
TMA effective interaction in the RMF Lagrangian, and describe pairing correlation 
by the density-independent delta-function interaction. We calculate single neutron 
and proton separation energies, quadrupole deformations, nuclear matter radii and 
differences between proton radii and neutron radii, and compare these results with 
the recent experimental data. 



1 Introduction 



One of the most fundamental quantity of the nucleus is its size, nuclear radius, 
which provides important information on the effect of the nuclear potential, 
shell effect and ground-state deformation. The presence of neutron skin in 
stable nuclei has been discussed since the late 1960s [1]. No evidence of appre- 
ciable neutron skin in stable nuclei has been observed. A clear evidence has 
not been obtained even for neutron-rich nuclei until 1990s. During the last 
decade, developments in accelerator technology and detection techniques [2] 
have led to an extended study of exotic nuclei near the limits of the (3- stabil- 
ity, especially those light nuclei with 4 < Z < 12. In particular, measurements 
of interaction cross sections (07) for light nuclei have been extensively studied 
at many radioactive nuclear beam facilities and provided important data on 
nuclear radii [3,4,6,5,7]. 
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In very neutron-rich nuclei, the weak binding of the outermost neutrons causes 
the formation of the neutron skin on the surface of a nucleus, and the forma- 
tion of one- and two-neutron halo structures. Measurements of Na isotopes 
[3,4] have revealed a monotonic increase in the neutron skin thickness with 
increasing neutron number. The established two-neutron halo nuclei are 6 He, 
n Li and 14 Be, and the one-neutron halo nuclei are n Be and 19 C. Recent data 
[5] present evidence for a one-neutron halo in 22 N, 23 and 23 F. On the proton- 
rich side evidence has been reported for the existence of a proton skin in 20 Mg 
[6]. More recently, Ozawa et al. [7] has confirmed for the first time a mono- 
tonic increase of proton skin. These observations indicate that neutron as 
well as proton skin are quite common phenomena in unstable nuclei far from 
the stability line. Data on ground-state deformations are also very important 
for the study of shell effect in exotic nuclei. In particular, they reflect major 
modifications in the shell structures, the disappearance of traditional and the 
occurrence of new magic numbers. The new magic number N — 16 has been 
reported in the light neutron drip line region [8]. Different ground-state de- 
formations of proton and neutron density distributions are expected in some 
nuclei with extreme isospin projection quantum number, T z . 

The Relativistic Mean Field (RMF) theory has been successfully applied to 
study global properties of nuclei from the proton drip line to the neutron 
drip line in the entire mass region [9,10]. The parameter sets, TM1 and TM2, 
were invented to describe ground-state properties of heavy nuclei (A > 40; 
TM1) and light nuclei (A < 40; TM2) [9]. The parameter set, TM1, has 
a favorable property of being able to reproduce the essential feature of the 
equation of state and the vector and the scalar self-energies of the relativistic 
Bruckner-Hartree-Fock theory for nuclear matter [11]. This parameter set was, 
however, not able to provide good quantitative results for light nuclei, and 
the introduction of the other parameter set, TM2, was forced to be made. 
Considering the RMF Lagrangian with the model parameters as completely 
a phenomenological model for description of nuclear ground states within the 
RMF approximation, Sugahara introduced a parameter set smoothly varying 
with mass and named it as TMA [12], which interpolates the TM1 and the 
TM2 parameter sets with very week mass dependence. We may interpret this 
mass dependence as a mean to effectively express the quantum fluctuations 
beyond the mean field level and/or the softness of the nuclear ground states 
in deformation, pairing and alpha clustering in light nuclei. On the other side, 
the parameter set, NL3 [13], is another very successful parameter set for the 
description of nuclear properties in the entire mass region, which originates 
from the original non-linear parameter set, NLI [14]. The NL3 parameter set 
has removed the unwanted features of NLI of not being able to obtain stable 
states for 12 C and 16 due to too large attractive contribution of the non-linear 
terms, while keeping the good properties of the NLI parameter set for heavy 
nuclei. It is, therefore, very interesting to compare the theoretical results of 
the RMF model with these parameter sets, TMA, TM2 and NL3 with the 
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newly obtained systematic data on the light mass region. 

In the present work, we apply the recently developed deformed RMF+BCS 
method with a density-independent delta-function interaction in the pairing 
channel [15] to the analysis of ground-state properties of Ne, Na, CI and Ar 
isotopes. As for the mean field part, we use the TMA, TM2 and the NL3 
parameter sets and compare, in particular, the proton and the neutron skins 
in Ar and Na isotopes. We compare also with the HF+BCS model with the 
Skyme parameter [16]. We note that there is a work by Lalazissis et al. [17] 
studying deformed light nuclei, Be, B, C, N, F, Ne and Na isotopes with the 
Relativistic Hartree-Bogoliubov (RHB) method. 

In Sect. 2, we shall present the RMF theory with deformation and pairing 
correlation. We provide the numerical results in Sect. 3 for Ne, Na, CI and Ar 
isotopes. Sect. 4 is devoted to the summery of the present study. 



2 Relativistic mean field theory with deformation and pairing 



Our RMF calculations have been carried out using the model Lagrangian 
density with nonlinear terms for both a and u mesons as described in detail 
in Ref. [15], which is given by 

C = - M)i/> + \dr<jfr(j - \mla 2 - \g 2 a* - \g^ - gjaifj 

-\R\ u R a ^ v + \m 2 pP lp a » - g p ^ llx r a ^(f a 



where the field tensors of the vector mesons and of the electromagnetic field 
take the following forms: 
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and other symbols have their usual meanings. Based on the single-particle 
spectrum calculated by the RMF method described above, we perform a state- 
dependent BCS calculation [18,19]. The gap equation has a standard form for 
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all the single particle states, i.e. 



A fc 



4e 



(3) 



where e k ' is the single particle energy and A is the Fermi energy. The particle 
number condition is given by 2 J] v\ = N. In the present work we use for the 



k>0 



pairing interaction a delta-function interaction, i.e. 
V = -V 5(r 1 -r 2 ), 



(4) 



with the same strength Vq for both protons and neutrons. The pairing matrix 
element for the delta-function interaction is given by 



V H = (n\V\jj) - (n\V\jj) = -V Jd 3 r UJ^J^.^j - 



(5) 



with nucleon wave function in the form 



Mr) 



\Wi(r) 



Xuit). 



(6) 



A detailed description of the deformed RMF+BCS method can be found in 
Ref. [15]. 

In the present study, nuclei with both even and odd number of protons (neu- 
trons) need to be calculated. We adopt a simple blocking method without 
breaking the time reversal symmetry. The ground state of an odd system is 
described by the wave function, 

a ] ki \BCS >= a f kl J[ (u k + v k a\a\) \vac> . (7) 



Here, \vac > denotes the vacuum state. The unpaired particle sits in the level 
ki and blocks this level. The Pauli principle prevents this level from partici- 
pating in the scattering process of nucleons caused by the pairing correlations. 
As described in Ref. [19], in the calculation of the gap, one level is "blocked": 

A k = -\ E , Vkk ' Ak ' • (8) 

2 k'^k^o ^ (e k , - \) 2 + A 2 k , 
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The level k\ has to be excluded from the sum because it can not contribute 
to the pairing energy. The corresponding chemical potential is determined by 



N = 1 + 2 




(9) 



fc^fci>0 



The blocking procedure is performed at each step of the self-consistent itera- 
tion. 



3 Ground state properties of light nuclei 

In the present work we take the mass-dependent parameter set TMA [9] for 
the RMF Lagrangian. Calculations with parameter set TM2 [9] and NL3 [13] 
are also performed for comparison. In the pairing channel, pairing strength 
V = 343.7 Mev fm 3 [15] is taken for both protons and neutrons in the 
density-independent delta-function interaction. For each nucleus, first the con- 
strained quadrupole calculation [20] is done to obtain all possible ground-state 
configurations; second we perform the non-constrained calculation using the 
quadrupole deformation parameter of the deepest minimum of the energy 
curve of each nucleus as the deformation parameter for our Harmonic Oscil- 
lator basis. In the case of several similar minima from the constrained cal- 
culation, we repeat the above procedure to botain the configuration with the 
lowest energy as our final result. The calculations for the present analysis have 
been carried out by an expansion in 14 oscillator shells for both fermion fields 
and boson fields. Convergence has been checked with more shells. Following 
Ref. [21], we fix Jtlu = 41 A^ 1 ^ 3 for fermions. In what follows we discuss the 
details of our calculations and the numerical results. 

3. 1 Single neutron and proton separation energy 

Single neutron separation energy: 



are very important and sensitive quantities of nuclei, where B(Z, N) are the 
binding energies of nuclei with proton number Z and neutron number N. We 



S n (Z,N) = B(Z,N)-B(Z,N-l) 



(10) 



and single proton separation energy: 
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plot the single neutron separation energies of our calculations with TMA set 
(empty square) for Ne, Na, CI and Ar isotopes, together with the available 
experimental values [22] (solid square) in Fig. 1. Good agreement between 
experiment and our calculations can be clearly seen. Calculations with TM2 set 
and NL3 set give practically the same results. For Ne isotopes, our calculations 
give relative small odd-even staggering than experiment except for 31 Ne. The 
most neutron-rich isotope ever observed for Ne isotopes is 34 Ne [23] and for 
Na isotopes it is 37 Na [23] , which are consistent with our present calculations 
(see also Table. 1 and Table. 2). The single proton separation energies of our 
calculations with TMA set (empty square) for Na and Ar isotopes, together 
with the available experimental values [22] (solid square) are plotted in Fig. 2. 
Once again, good agreement between experiment and our calculations can be 
clearly seen. For 21 Na, the difference between experiment and our calculations 
seems somewhat large. We should note that, from the definition of S p (Eq. 11), 
the single proton separation energy for 21 Na is related to the binding energies 
of 21 Na and 20 Ne. While in our calculations the calculated binding energy for 
20 Ne is about 1.3 MeV smaller than experiment and for 21 Na is about 1.2 
MeV larger than experiment. Thus a total difference of 2.5 MeV is given to 
the single proton separation energy for 21 Na. For 35 Ar, the relatively large 
difference comes from the fact that the calculated binding energy for 34 CI 
is about 2.5 MeV larger than experiment, while calculated binding energy 
for 35 Ar is about 1 MeV larger than experiment. The last bound proton-rich 
isotope in Na isotopes is predicted to be 20 Na in our calculations, which has 
a single proton separation energy 3.099 MeV. Its one neutron less neighbor, 
19 Na, has almost a zero single proton separation energy -0.013 MeV (see also 
Table. 1 and Table. 2). The last bound proton-rich isotope in Ar isotopes is 
predicted to be 31 Ar in our calculations, whose S p is 0.741 MeV (see also 
Table. 3 and Table. 4). All these predictions are consistent with available 
experimental knowledge [24]. 

In conclusion, our calculations give a very good description of single neutron 
and proton separation energies for Ne, Na, CI and Ar isotopes. All three cal- 
culations with different parameter sets, TMA, TM2 and NL3, give essentially 
the same results. Therefore, we can conclude that the binding energies and sin- 
gle neutron (proton) separation energies are reproduced very well with TMA 
and NL3 parameter sets, which have been checked by many observables in the 
entire mass region. 



3.2 Root mean square nuclear radius 

The root mean square neutron, proton, and matter radii are other important 
basic physical quantities for nuclei in addition to the single neutron (proton) 
separation energies. In the RMF theory, the root mean square (rms) neutron, 
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proton and matter radii can be directly deduced from the neutron, proton and 
matter density distributions, p n , p p and p m , 

Ri = (r?)1/2 = { /^ }V2 _ (12) 



where the index % (= n,p,m) denotes the corresponding neutron, proton and 
matter density distributions. First let us have a look at the rms neutron and 
rms matter radii of Na and Ar isotopes. In Fig. 3, we plot results of our calcu- 
lations with TMA set and the available experimental values for Na [3,4] and 
Ar isotopes [7]. For neutron radii and matter radii of Na isotopes, we notice 
that our calculations agree very well with experiment except for 22 Na, which 
has recently been attributed to an admixture of the isomeric state in the beam 
[4]. For Ar isotopes, although the mass dependence of the experimental data 
are reproduced quite well, some discrepancies remain. First, the theoretical 
predictions are at the upper limits of experimental values. Calculations with 
TM2 set and NL3 set are not better either. Second, the larger radii of Z = N 
nuclei 36 Ar than those of its neighbors has been interpreted to be a possi- 
ble alpha-cluster structure [25]. Third, our calculations do not reproduce the 
abrupt decrease of 37 Ar and 38 Ar. The relatively larger radii of proton drip line 
nucleus 31 Ar can be understood easily if we note that 31 Ar is the last bound 
proton-rich isotope in Ar isotopes. Second we plot the charge isotope shifts of 
our calculations with TMA set (empty square) together with the experimental 
values (solid square) for Na isotopes [26] and Ar isotopes [27] in Fig. 4. Cal- 
culations with TM2 set (solid circle) and results of HF+BCS model [16] are 
also shown for comparison. For Na isotopes, significant difference between the 
theoretical predictions and the experimental values exist in two regions. The 
first is 25 Na, and the second is around the region of neutron drip line A > 29. 
While for Ar isotopes, calculations with TM2 set agree very well with the 
experimental values within the error bar except for 40 Ar. For both isotopes, 
results with TM2 set are closer to the experimental values. It is not surprising 
because TM2 is a parameter set specially made for light nuclei. For further 
comparison with other theories, we include the results of the HF+BCS model 
[16] in Fig. 4. As for Ne isotopes, the isotope shifts start to deviate from 26 Ne 
to the lower side except for 31 Ne. On the other hand, the results for Ar isotopes 
agree very well with experiment and are similar to those of the RMF+BCS 
with TM2 parameter set except for 32 Ar. 

To conclude, the deformed RMF+BCS method describe the rms neutron radii 
and matter radii very well. While for rms charge (proton) radii, although 
the basic trend is reproduced quite well, the results are not so satisfactory 
considering that we obtain the parameter set by fitting the charge radii of 
certain nuclei. However, we might have to pay more attention to the various 
cluster phenomena, which make it difficult to apply the mean field theory to 
light nuclei. 



7 



3. 3 Neutron skin and proton skin 



Different theories, both relativistic and non-relativistic ones, have predicted 
the existence of proton skin and neutron skin since long time ago, but only 
recently experimental physicists have proved the theoretical predictions. In 
Fig. 5, the proton skin AR = R p — R n is plotted against the difference between 
the proton and neutron separation energy AS = S p — S n . It can be clearly 
seen that AR has a strong correlation with AS*. Thus, it is shown that the 
difference between the proton and neutron Fermi energy is the driving force 
for the creation of the skin phenomenon. Such a correlation has already been 
predicted by RMF theory since a decade ago [28]. For Na isotopes, calculations 
with NL3 set predict smaller AS for 28-31 Na, while calculations with TMA 
and TM2 sets are in better agreement with experimental AS values. For Ar 
isotopes, all three calculations failed to reproduce the experimental observed 
proton-skin in 37 Ar and 38 Ar. In Fig. 6, we plot neutron skin R n — R p against 
mass number A for Ne, Na, CI and Ar isotopes. From Figure. 6, it is quite clear 
that neutron (proton) skin are quite common for neutron-rich (proton-rich) 
nuclei. 



In the Relativistic Mean Field theory, the formation of skin or halo is usually 
explained by the picture of a core plus a few occupied loosely bound (unbound) 
states of small angular momentum. The halo in n Li has been successfully re- 
produced in this picture [29] where the occupations of lpi/2 and 2si/ 2 state are 
found to be important. Based on the similar picture, giant-halos in 124-138 Zr 
have been predicted by both RCHB [30] method and resonant RMF-rBCS 
method [31]. In Ref. [31], it has been shown that 3p3/2, 2f 7 / 2 and 3pi/ 2 states 
contribute most to the formation of halo in 124 Zr. In order to reproduce the 
sudden increase of rms neutron (proton) radii, which is one indication of halo 
structure, the Relativistic Mean Field model must be solved in coordinate 
space. Although, compared with RCHB method, quite similar two neutron 
separation energies are obtained in the deformed RMF+BCS method [15], the 
sudden increase of rms neutron radii for 124 ~ 138 Zr are not reproduced well. 
Some improvements on the expansion method could improve the tail behavior 
of resonant wave functions and solve this problem finally. Possible candidates 
are the expansion method in the local transformed Harmonic Oscillator basis 
[32] and the expansion method in a Woods-Saxon basis [33] . This work is under 
way now. However, as we can see from our calculations, such a disadvantage 
of Harmonic Oscillator basis does not influence our conclusions here. 
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3.4 Deformation parameter 

Compared with binding energies and nuclear radii, predicted deformation pa- 
rameters of nuclei are quite different from model to model. One reason is, of 
course, the possible shape coexistence. Another reason is that deformation pa- 
rameters are more sensitive to model details. On the experimental side, we can 
obtain the information of nuclear deformation from measurements of B(E2) f 
values. The B(E2) j values are basic experimental quantities that do not de- 
pend on nuclear models. Assuming a uniform charge distribution out to the 
distance R{0,4>) and zero charge beyond, (3 2p [34,15] is related to B(E2) j by 

(3 2p = (An/3ZR 2 )[B(E2) | /e 2 ] 1 / 2 , (13) 

where Rq has been taken to be 1.2A 1 / 3 fm and B(E2) | is in units of e b . 
B{E2) | values have been measured mainly for some even-even nuclei. 

In Fig. 7, we plot the deformation parameter f3 2p i m ) f° r Ne, Na, CI and Ar 
isotopes against mass number A. For our RMF+BCS calculations, we plot (3 2p 
because as we can see from Eq. (13) that what we obtain from B(E2) j is 
actually /3 2p , not (3 2m . For the Finite Range Droplet Model (FRDM) [35], the 
plotted deformation parameter is (3 2m . Except for exotic nuclei with extreme 
N/Z ratios, we note that these two values are very close to each other (see 
also Table. 1 ~ 4). From Fig. 7, we see that experimental values show that Ne 
isotopes are super deformed {(5 2p ~ 0.6) in the region A = 18 ~ 28. Except for 
20 Ne and 22 Ne, our calculations failed to reproduce 18 Ne, 24 Ne, 26 Ne and 28 Ne. 
For Na isotopes, both our calculations and FRDM predict prolate shapes. For 
CI isotopes, a transition between prolate shapes and oblate shapes have been 
predicted by both methods. For Ar isotopes, our calculations predict oblate 
shapes for most nuclei while FRDM predict spherical shapes. Except for 38 Ar, 
our predictions agree very well with experimental values. We should mention 
here that, because from the B(E2) j values, we can not distinguish prolate 
or oblate shapes, we just put a random sign before the experimental value 
to make it close to our predictions. This does not change the essence of our 
comparison. 



4 Summary 

We have studied light nuclei, in particular, Ne, Na, CI and Ar isotopes, in the 
framework of the deformed RMF+BCS method. The deformation is treated 
by using the expansion method in the deformed Harmonic Oscillator basis. We 
have used TMA, TM2 and NL3 effective interactions in the RMF Lagrangian. 
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In addition we have treated the pairing correlations in terms of a density- 
independent delta-function interaction. 

We have calculated neutron (proton) separation energies, quadrupole deforma- 
tions, nuclear neutron (proton) and matter radii. All the results are presented 
in the form of tables (Table. I ~ IV). We have compared the calculated results 
with three parameter sets, TMA, TM2 and NL3, with the available experi- 
mental values. Both TMA and NL3 parameter sets are often used in the RMF 
theory in the entire mass region from the proton drip line to the neutron drip 
line. We have used TM2 parameter set also, which is extracted particularly for 
light nuclei. All the three parameter sets provide similar results in our present 
calculations. 

The single neutron (proton) separation energies are reproduced very well 
within the deformed RMF+BCS method with all the three parameter sets. 
The general trends of the rms neutron radii come out to be also very good. 
However, the rms proton radii seem to be generally overestimated in the de- 
formed RMF+BCS method in comparison with experimental values. We have 
extracted also the nuclear deformations by the minimization method. In gen- 
eral, the calculated results are in good agreement with experimental values 
extracted from the BE (2) j values. The discrepancy on the proton radii for 
these light nuclei may be related with the softness of the light nuclei in de- 
formation and in alpha-clustering, which are not included in the mean field 
Lagrangian. We have to study the effect of the softness in future works 
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Fig. 1. The single neutron separation energis, S n , of Ne, Na, CI and Ar isotopes. 
Results obtained from the deformed RMF+BCS calculations with TMA set are 
compared with available experimental data [22]. 
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Fig. 2. The single proton separation energies, S p , of Na and Ar isotopes. Results 
obtained from the deformed RMF+BCS calculations with TMA set are compared 
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Fig. 3. The rms neutron radii, R n , and rms matter radii, R m , of Na and Ar isotopes. 
Results obtained from the deformed RMF+BCS calculations with TMA set are 
compared with available experimental data [3,4,7]. 
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Fig. 4. The charge isotope shifts of Na and Ar isotopes. Results obtained from the 
deformed RMF+BCS calculations with TMA set and TM2 set are compared with 
available experimental data [26,27] and results of HF+BCS model [16]. 
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with all the three parameter sets, TMA, TM2 and NL3, from top to bottom. The 
corresponding experimental data are taken from Ref. [3,4,7,22]. 
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Fig. 6. The neutron skins, R n — R p , of Ne, Na, CI and Ar isotopes. Results ob- 
tained from the deformed RMF+BCS calculations with TMA set are compared 
with available experimental data [3,4,7]. 
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Table 1 

The ground state properties of Ne isotopes (Z = 10) calculated with the parameter 
set TMA. Listed are the total binding energy, B, the binding energy per nucleon, 
B/A , charge, neutron, proton, and matter root mean square radii, R c , R n , R p 
and R m , and the quadrupole deformation parameter for the neutron, proton and 
matter distributions, fa n , fap an d fa m , with A the mass number and N the neutron 
number. 
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Table 2 



The same as Table 1, but for Na isotopes (Z = 11). 
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Table 3 

The same as Table 1, but for CI isotopes (Z = 17). 
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Table 4 



The same as Table 1, but for Ar isotopes (Z = 18). 
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